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ABSTRACT 

There is increasing evidence that supermassive black holes in active galactic nuclei 
(AGN) are scaled-up versions of Galactic black holes. We show that the amplitude of high- 
frequency X-ray variability in the hard spectral state is inversely proportional to the black 
hole mass over eight orders of magnitude. We have analyzed all available hard-state data from 
RXTE of seven Galactic black holes. Their power density spectra change dramatically from 
observation to observation, except for the high-frequency (> 10 Hz) tail, which seems to have 
a universal shape, roughly represented by a power law of index -2. The amplitude of the tail, 
Cm (extrapolated to 1 Hz), remains approximately constant for a given source, regardless of 
the luminosity, unlike the break or QPO frequencies, which are usually strongly correlated 
with luminosity. Comparison with a moderate-luminosity sample of AGN shows that the am- 
plitude of the tail is a simple function of black hole mass, Cm — C/M, where C rs 1.25 M Q 
Hz -1 . This makes Cm a robust estimator of the black hole mass which is easy to apply to low- 
to moderate-luminosity supermassive black holes. The high-frequency tail with its universal 
shape is an invariant feature of a black hole and, possibly, an imprint of the last stable orbit. 
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1 INTRODUCTION 

Astrophysical black holes are very simple objects, completely char- 
acterized by their mass and spin. Hence, the gravitational poten- 
tial around a black hole simply scales with its mass. An important 
question in high-energy astrophysics is whether the accretion flow 
properties scale with the black hole mass in a simple manner, or, 
more specifically, whether active galactic nuclei (AGN) are scaled- 
up versions of Galactic black hole binaries (BHB). 

One of the ways of tackling this problem is to study X-ray 
variability, which is observed in accreting black holes of all masses. 
Recent advances in mass estimates of AGN central black holes lead 
to discovery of dependence of the observed variability properties 
on mass. Long X-ray monitoring campaigns allowed to construct 
power density spectra (PDS) of accreting supermassive black holes 
which turned out to have a roughly of (broken) power-law shape. 
The variability amplitude (the excess variance; e.g. Lu & Yu 2001; 
Markowitz & Edelson 2001) and the frequency of the break (e.g. 
M c Hardy et al. 2004, 2006) can depend on the black hole mass. 
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In order to use the X-ray variability for mass measurement 
we need a property which scales only with the black hole mass, 
and does not change with accretion rate. The break frequency 
does not satisfy this condition, as it changes significantly with the 
accretion rate, in X-ray binaries (e.g. Done & Gierlinski 2005). 
M c Hardy et al. (2006) showed that a more general relation holds 
between the break frequency, vi, and the black hole mass, M: 
Vh — AL^ ol /M c ', where A, B and C are constants. This relation 
includes a significant dependence on the source bolometric lumi- 
nosity, Lboi- 

It was already suggested by Hayashida et al. (1998) that mea- 
suring the normalization of the high-frequency tail of the power 
spectrum, well above the high-frequency break, is an interesting 
possibility for black hole mass measurement. Equivalently, one 
can use the excess variance, cr^xs* measured for short data sets. 
This general line was followed by Czerny et al. (2001), Papadakis 
(2004), Nikolajuk, Papadakis & Czerny (2004, hereafter N04) and 
Nikolajuk et al. (2006, hereafter N06). However, the method was 
not reliably checked against the dependence on the source accre- 
tion rate. 

In this paper we put the idea of <j%xs ^ M~ x correlation to 
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the test. We use an extensive set of BHB observations to see if and 
when on XS is constant for a given mass and whether it anticorre- 
lates with the black hole mass. 



2 HIGH-FREQUENCY POWER 

Power density spectra of many AGN can be approximated by a bro- 
ken power law, with power P v oc u^ 1 below and P v oc u~ 2 above 
the break frequency, Vb (e.g. Markowitz et al. 2003b), where P v 
is the power spectral density normalized to the mean and squared. 
A second break at lower frequencies, below which the power is 
roughly P v oc v°, has been also observed (e.g. Pounds et al. 2001; 
Markowitz, Edelson & Vauhgan 2003a). At the zeroth order of ap- 
proximation this is consistent with the PDS observed in stellar- 
mass BHB in the hard spectral state. Fig. Q] shows a sample of 
PDS from Galactic BHB in the hard state (details of the data re- 
duction are described in Section[3]l. Plainly, these spectra are much 
more complex that a doubly-broken power law, with multiple broad 
and narrow noise components, usually well described by a series of 
Lorentzians (e.g. Pottschmidt et al. 2003). However, despite this 
complexity, the entire spectral shape roughly resembles a (doubly) 
broken power law. 

N04 assumed that the break frequency is inversely propor- 
tional to the black hole mass, while the P v oc v~ part of the 
PDS below the break (the 'flat top' in vP v diagrams) has constant 
normalization, independent of the black hole mass. Yet inspection 
of several BHB power spectra clearly shows that neither of these 
is constant for a given source. The break frequency is known to 
change with accretion rate (e.g. Done & Gierlihski 2005). The 'flat 
top' normalization can change as well, as one can see in GX 339^1 
spectra in Fig.Q] 

There is, however, one feature of these power spectra that 
remains remarkably invariant: the high-frequency spectral shape, 
above Vb- For a given source it can be roughly described by a single 
power law with constant index of 1.5-2.0, and constant normaliza- 
tion for various observations differing in luminosity by more than 
one order of magnitude. In this paper we test the idea that the high- 
frequency part of the PDS remains fairly constant for a given source 
and scales with the black hole mass. This is a simple refinement of 
the idea proposed by N04 and later developed by N06. 

Here we do not make any assumptions about how the char- 
acteristic frequencies (e.g. break frequency) depend on black hole 
mass. Instead, we assume that the the PDS above the break fre- 
quency (the high-frequency tail) has a universal spectral shape 
(roughly oc v~ 2 ) with normalization depending on the black hole 
mass. This can be written as 

P v = C M {v/v )- 2 , (1) 

where va is an arbitrary frequency which we chose to be z^o = 1 Hz. 
Thus, Cm (in units of Hz - ) is the normalization of the (extrapo- 
lated) high-frequency tail at 1 Hz. 

The assumption that Cm is unique function of the black hole 
mass would directly correspond to the original assumption of N04 
about constancy of P(yb)Vb if Vb were constant for a given black 
hole mass. Due to limited statistics it is often difficult to study de- 
tails of the high-frequency shape of the PDS. Therefore, we sim- 
plify the situation by calculating the amplitude, or the excess vari- 
ance, of variability in a given frequency band significantly above 
the break. This can be done directly from a light curve or by inte- 
grating the PDS. The excess variance calculated between frequen- 
cies v\ and V2 (both greater than Vb) is: 



Source Name 


Start 


End 


XTEJl 118+480 


2000-03-29 


2000-08-08 




2005-01-13 


2005-02-26 


4U 1543^17 


2002-06-17 


2002-10-11 


XTEJ1550-564 


1998-09-07 


1999-05-20 




lfi(i(| A/1 1 A 

2UUU-U4-1U 


ZUUU-U /-to 




2001-01-28 


2001-04-29 




2002-01-10 


2002-03-05 




2003-03-27 


2003-05-16 


XTE J1650-500 


2001-09-06 


2002-04-21 


GRO J1655-40 


2005-02-20 


2005-11-11 


GX 339^1 


2002-04-02 


2003-05-06 




2003-12-28 


2005-08-12 


Cyg X-l 


1996-02-12 


2006-01-12 



Table 1. Log of RXTE observations. Each set of data corresponds to one 
transient outburst. For Cyg X- 1 we used all data publicly available in Febru- 
ary 2007. 



0"nxs = / Pvdv = CmVo { — ) ■ (2) 

/ V ^1 V2 J 

The key assumption, which we want to test in this paper, is that 
Cm is inversely proportional to the black hole mass, Cm — C/M, 
where C is a constant. We note that our C is the same constant as 
constant C defined by N04 in eq. 4, divided by v 2 . 



3 DATA REDUCTION AND SELECTION 

We used publicly available Rossi X-ray Timing Explorer (RXTE) 
Proportional Counter Array (PCA) data of seven black hole bi- 
naries, listed in Table Q] First, we extracted background-corrected 
energy spectra from Standard-2 data (top layer, detector 2 only) 
for each pointed observation. These were used to create hardness- 
intensity diagrams (HID). Intensity is defined as the total 2-60 keV 
count rate and the hardness ratio is the ratio of count rates in energy 
bands 6.3-10.5 and 3.8-6.3 keV. We also calculated power density 
spectra (PDS) for each observation from full-band (2-60 keV) data 
in 0.0039-128 Hz frequency band. We subtracted the Poissonian 
noise from the PDS, corrected them for dead- time effects (Revnivt- 
sev, Gilfanov & Churazov 2000) and background (Berger & van 
der Klis 1994) 

We used HIDs to identify hard state spectra. In transients these 
spectra lie on the vertical branch of the diagram at hardness ra- 
tio >0.8, corresponding to the rise or decay of the outburst. While 
in GX 339—4 the transition from the vertical hard branch into in- 
termediate horizontal branch is abrupt and well defined, in XTE 
J1550-564, XTE J1650-500 and GRO J1655-40 we also included 
a few points from the part of the HID where the hard branch turns 
into intermediate, as they have sufficiently large hardness ratio to 
be (potentially) classified as the hard state. In Cyg X-l we arbitrar- 
ily selected observations with hardness ratio greater than 0.9. We 
also arbitrarily rejected all data with average count rate less than 
10 counts per s per PCU as they did not have enough statistics to 
robustly calculate the high-frequency power. All the selections are 
shown in Figs. [2] and [3] 

For a given source (and separately for each outburst in case of 
GX 339-4 and XTE 11550-564) we measured the high-frequency 
power, Onxs> with their statistical errors, a, for each pointed obser- 
vation in the hard state. The power (or excess variance) was mea- 
sured by integrating the PDS over the frequency band 10-128 Hz 
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Figure 1. Power density spectra of black hole X-ray binaries. Each panel shows three spectra in the low/hard spectral state, selected from the data analyzed 
in Sec. [3] The selection demonstrates that despite dramatic changes in the PDS shape at lower frequencies (including variable QPO), the high-frequency part 
( >10 Hz) of the PDS remains relatively constant. 



(eq.[2j. Then, we computed the mean amplitude, (Cm) for a given 
source (or a particular outburst of the source) weighted by errors 
of Cm- The error of (Cm) was estimated from \ 2 statistics for 
A\ 2 = 2.7, i.e. corresponding to 90 per cent confidence limits. 

Despite the initial count rate selection, some data sets (in par- 
ticular from short observations) gave high-frequency power with 
large errors. We discarded all these, setting an arbitrary upper limit 
on error, a < O.Scr^xs- 



4 RESULTS 

Figs.[2]and[3]show selection of hard-state data (hardness-intensity 
diagrams on the left) and measured amplitude of the high- 
frequency tail, Cm in the right-hand panels. It is obvious from these 
diagrams that Cm is not constant for a given source and varies from 
one pointed observation to another. These variations are not very 
significant, though. With very few exceptions (2 observations of 
GRO J1655^t0 and 3 of GX 339^1) individual Cm measurements 
are within 3<r of the mean. The dispersion is higher in Cyg X-l 
where statistics is better and errors on Cm smaller. Generally, we 
find that Cm does not significantly depend on the source brightness 
(count rate). 

More significant differences can be found for different out- 
bursts of a given source. We have analyzed five outbursts of XTE 
J1550-564. The three hard-state outbursts between 2001 and 2003 
(see Table [2j were similar in all properties, so we analyzed them 
together. They yield mean (Cm) = 0.101 ± 0.002 Hz" 1 . The hard 
state in 2000 outburst gave higher (Cm) = 0.137 ± 0.002 Hz" 1 . 
However, the 1998 outburst was very different, with much lower 
and quickly changing high-frequency power, Cm ~ 0.05 Hz -1 . 
We have excluded the onset of 1998 outburst from further analysis, 
as it might have represented a different accretion state (we discuss 



Source Name Outburst (Cm) (Hz 1 ) Thard (days) 



XTE Jill 8+480 


2000 


0.098 ± 0.005 


130* 




2005 


0.14 ±0.02 


26* 


4U 1543^17 


2002 


0.14 ±0.01 




XTEJ1550-564 


1998 


~0.05 


3 




2000 


0.137 ±0.002 


14 




2001 


0.095 ± 0.001 


90* 




2002 


0.104 ±0.004 


50* 




2003 


0.101 ± 0.002 


50* 


XTEJ1650-500 


2001 


0.113 ± 0.001 


8 


GRO J 1655^10 


2005 


0.179 ± 0.003 


12 


GX 339^1 


2002 


0.088 ± 0.001 


20 




2004 


0.121 ± 0.003 


180 


Cyg X-l (hard) 




0.0868 ± 0.0003 




Cyg X-l (soft) 




0.139 ± 0.004 





Table 2. Mean amplitude of the high-frequency power, (Cm) (eq[D- We 
show data for each outburst separately. Thard is approximate duration of 
the initial hard state during the rise of the outburst. Asterisks denote hard- 
state only outbursts. 4U 1543^47 was not observed in the initial hard state 
and Cyg X-l is a persistent source. 



this in detail in Sec.0. Similarly, GX 339-4 and XTE Jl 1 18+480 
showed different Cm during different outbursts, though at least in 
the case of the latter one this could have been due to systematic 
effects at low count rates, which we discuss later in this section. 

We looked into the dependence of Cm on the hardness ratio, 
which can be regarded as a crude indicator of the spectral state. We 
did not find any clear general trend, as illustrated in Fig. [4] 

Additionally, we found variability amplitude for three other 
black hole candidates with no black hole mass estimates. The rel- 
atively high value of (Cm) = 0.15t^ Hz" 1 obtained for XTE 
J1720-3 18 seems to be consistent with a rather low black hole mass 
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Figure 2. Panels on the left show hardness-intensity diagrams. Solid line 
shows the entire outburst of a given source, while the data points correspond 
to the hard-state selection of pointed observations, used for further analy- 
sis. Diagrams on the right show the high-frequency variability amplitude, 
Cm (eq. E}> as a function of count rate, for selected hard-state observa- 
tions. The horizontal line shows the weighted average of Cm f° r the entire 
outburst. XTE Jl 118+480 shows data from two outbursts: 2000 (red trian- 
gles) and 2005 (black filled squares). XTE J1550-564 diagram shows data 
from five outbursts: 1998/1999 (red open triangles), 2000 (blue squares) 
and 2001, 2002 and 2003 together (green crosses). The circled point in XTE 
Jl 118+480 panel was used to estimate background effects (see Section[4). 



of 5 Mq, as also suggested by Cadolle Bel (2004) from disc spec- 
tral fitting. For H1743-322 (=XTE J1746-322) and XTE J1748- 
288 we found (C M ) = 0.090 ± 0.015 and 0.056 ± 0.008 Hz -1 , 
respectively. 

The hard X-ray spectral state, which we study in this paper, 
is dominated by Comptonized emission. It would be very inter- 
esting to see if similar behaviour can be seen in the other spec- 
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Figure 3. Continuation of Fig. [2] GX 339^4 diagrams shows two outbursts: 
2002/2003 (black squares) and 2004 (red triangles). 



tral state dominated by Comptonization, the very high (or steep 
power law) state. We have selected the very high state data from 
hardness-intensity diagrams of XTE J1550-564 and GRO J1655- 
40, as shown in Fig. [5] These observations have significantly less 
variability at frequencies > 10 Hz then the hard state data, so the re- 
sulting Cm is small (black squares in Fig. [5}. But we must remem- 
ber that the soft X-ray part of the spectrum (where PCA has high- 
est sensitivity) has strong contribution from the cold accretion disc, 
which can dilute the variability coming from Comptonization (e.g. 
Done & Gierlinski 2005) hence suppressing the observed power. 
Therefore, we extracted additional power spectra at higher ener- 
gies (above PCA channel 36, roughly corresponding to energy of 
14 keV), where the disc influence is negligible. The high-frequency 
variability amplitude from these data is much higher, as shown in 
blue triangles in Fig. [5] Though there is a large scatter in Cm from 
individual observations, the mean {Cm) is of order of 0.1 Hz -1 in 
both sources, consistent with the hard state results. Therefore, the 
high-frequency variability from Comptonization appears to be very 
similar in both hard and very high states (but see discussion in Sec. 
0. 

The soft-state data from our BHB sample are strongly domi- 
nated by the disc emission. The only source with reasonable count 
rates at higher energies in the soft state is Cyg X-l. We applied 
the same approach to Cyg X-l soft-state data (arbitrarily chosen 
hardness ratio in the range 0.45-0.55), as to the very high state 
data in the previous paragraph. The result can be seen in Fig. 
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Figure 4. Dependence of amplitude of high-frequency variability, Cm > on 
the hardness ratio. There is no clear trend of Cm as a function of spectral 
states. 



[6] The scatter in individual data points in significant, with mean 
(Cm) = 0.138 ± 0.004 Hz^ 1 . Since the fit of the constant to the 
data is rather poor, xi = 156/73, the error on (Cm) is indicative 
only. 
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Figure 5. Same as in Figs.|2]and[3] but for the very high (steep power law) 
spectral state. Black filled squares represent Cm calculated from 2-60 keV 
PCA data, grey (blue in colour) open triangles show the high-energy data 
(roughly 14-60 keV). 
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Figure 6. Same as in Fig. [3] but for the soft spectral state of Cyg X-l. 



5 CAVEATS 

We tested our results for possible systematic effects. Many of our 
observations have low count rate, which can affect the resulting 
power spectra. At low count rates, PCA channels 0-7 can cre- 
ate artificial power due to problems with the on-board computer 
(Revnivtsev, private communication). To test the effect of this we 
have removed PCA channels 0-7, where PCA configuration al- 
lowed for that, and calculated new values of Cm- We found that 
this had a negligible effect on our results. 

The white noise level was subtracted from the power spectra 
during data reduction. If the white noise level was not estimated 
correctly, this could have influenced the amplitude of variability 
and the resulting Cm- To test this we have selected observations 
where high timing resolution was available and calculated power 
spectra up to 1024 Hz. As we do not expect any significant power 
above a few hundred Hz from black holes (Sunyaev & Revnivtsev 
2000), we assumed that the 512-1024 Hz power could be used as 
a good white noise estimator. We reanalyzed these data and calcu- 



lated new values of Cm- Again, the effect on our results turned out 
to be negligible. 

Background effects can be potentially important for estimat- 
ing the amplitude of variability, which is defined as (rms/mean) 2 . 
Our power spectra were calculated from light curves not corrected 
for background, so their power is [rms/(R a + Rb)] 2 , where R s 
and Rt are source and background mean count rates, respectively. 
These PDS were then multiplied by [(R s +Rb) / Rs] 2 , where source 
and background count rates were estimated from light curves. In 
the dimmest observations R s and Rb are comparable. To estimate 
possible background inaccuracy (which is modelled in the PCA 
rather then measured) we extracted one PCA and HEXTE spec- 
trum (using standard HEASARC reduction techniques, adding 1 
per cent systematic errors in the PCA) of XTE Jl 118+480, cor- 
responding to a high-CM and low count rate point circled in Fig. 
[2] (ID 90111-01-02-07, observed on 2005-01-24). We have fitted 
the joined PCA/HEXTE spectrum with a simple Comptonization 
model (see Done & Gierliriski 2003 for details), using PCA in 
3-40 and HEXTE in 20-200 keV band. The fit was good ( X 2 = 
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Source Name 


M(M ) 


Reference 


XTEJ1118+480 


8.5 (7.9-9.1) 


Gelino et al. (2006) 


4U 1543^7 


9.4 (7.4-11.4) 


Park et al. (2004) 


XTEJ1550-564 


10(9.7-11.6) 


Orosz et al. (2002) 


XTEJ1650-500 


5 (2.7-7.3) 


Orosz et al. (2004) 


GRO J1655-40 


6.3 (5.8-6.8) 


Greene, Bailyn & Orosz (2001) 


GX 339^* 


6 (2.5-10) 


Cowley et al. (2002) 


Cyg X-l 


20 (13.5-29) 


Ziolkowski (2005) 



Table 3. Black hole masses in X-ray binaries, used in this paper. 



137/145) with no strong residuals. Then we changed the level of 
PC A background by ± 10 per cent. The fit with 90 per cent of back- 
ground was only marginally worse (X 2 = 142/145), while the 110 
per cent background resulted in a rather poor fit {\ 2 = 171/145). 
In both cases there were strong residuals in the PCA and disagree- 
ment with HEXTE data above -~25 keV. This shows that the back- 
ground in the low count rate data is estimated with accuracy much 
better then 10 per cent. Hence, the uncertainty on Cm due to back- 
ground estimation is no more than a few per cent. The increase of 
Cm by factor 2-5 at low count rates seen in Fig.|2]cannot be caused 
by incorrect background. Since we see this effect in most sources 
below the same count rate of ~20 s" 1 (regardless of the distance, 
hence at different luminosities) it must be of (unknown) instrumen- 
tal origin. We would like to stress, however, that the increase is not 
statistically significant, typically less then 2a. 

The high-frequency variability is known to depend on energy 
in some sources (e.g. Nowak et al. 1999). This is important when 
comparing BHB with AGN, as we look at different parts of the 
Comptonized spectrum, AGN data showing higher scattering or- 
ders than BHB. We have tested our data for energy dependence. 
This was possible only in bright observations from XTE J1550- 
564, GX 339-4 and Cyg X-l, where statistics at higher energies 
was good enough. The high-energy (above ~14 keV) data give 
(Cm) = 0.13 ± 0.01 Hz" 1 for the 2000 outburst of XTE J1550- 
564 and (Cm) = 0.09 ± 0.01 Hz" 1 for the 2002 outburst of GX 
339-4. These values are consistent with the broad-band data, which 
in the PCA is dominated by soft X-rays (see TableO, which sug- 
gests that the high-frequency amplitude is not energy-dependent in 
these sources. On the other hand, similar approach to Cyg X-l gave 
(Cm) = 0.109 ± 0.001 Hz" 1 , higher by about 25 per cent higher 
then the broad-band data. 



6 DEPENDENCE ON MASS 

Fig. 0a) shows the dependence of Cm on black holes mass. We 
used best currently available mass estimates of black holes in X-ray 
binaries, as summarized in Table[3] XTE J1650-500 does not have 
good mass estimate, though Orosz et al. (2004) found an upper limit 
of 7.4 M©. We assumed the mass function, 2.7 Mo, as the lower 
limit and we adopted the actual mass in the middle of this interval, 
at 5 M . We show data from different outbursts of XTE 1 1 1 8+480, 
XTE J1550-564 and GX 339-4 separately. 

Clearly, there is no apparent correlation between black hole 
mass and Cm , though we must bear in mind that mass estimates in 
X-ray binaries are not very accurate. Moreover, different outbursts 
giving slightly different (Cm) create additional dispersion. Hence, 
the overall uncertainties are rather large, so the expected relation 
Cm = C/M cannot be robustly confirmed from X-ray binaries. 
To do this, we need to extend our studies to supermassive black 



holes. N06 compared masses of a sample of Seyfert 1 galaxies mea- 
sured by reverberation method with masses from high-frequency 
variability, using the value of constant C derived from Cyg X-l 
observations. Here we use the same sample of AGN in order to 
compare them with our much larger sample of BHB, constrain the 
Cm = C/M relation better and get a 'big picture' overview of 
variability properties for all masses of black holes. In Fig.[7Jc) we 
plotted the sample of of N06; panel b shows the overview of stellar- 
mass and supermassive black holes. 

The red diagonal line in the diagrams represents the best- 
fitting function Cm = C/M, with C = 1.25±0.06 M Q Hz" 1 . 
We would like to point out that this particular value depends on the 
selection of X-ray binary data, as different outbursts can give dif- 
ferent (Cm). Also, errors on mass are non-Gaussian in many cases, 
so the error on C, given here, is indicative only. It is interesting to 
notice that the constant C is 1.24±0.06 M Q Hz -1 for X-ray bina- 
ries only and 1.5lt°jg Mq Hz 1 for AGN only; both values are 
consistent within error limits. We also plotted a line (green dashed) 
corresponding to the soft state of Cyg X-l, for comparison (with 
C = 2.77 M Hz" 1 ). Some of the AGN are consistent with the 
soft-state line. 

We also tested whether the relation between variability and 
mass is really linear. We fitted a more general form of a power-law 
dependence, Cm = C/fM/Mg)"" to all BHB and AGN data and 
found C f = 1.19 ± 0.05 Hz" 1 and index a = 0.98 ± 0.01 very 
close to linear relation. Hence, we opt for the simpler solution and 
regard the linear relation as well established. 

The Cm = C /M relation is in excellent agreement with the 
data spanning over eight orders of magnitude in mass. It gives 
us robust confirmation of our hypothesis that the high-frequency 
power is inversely correlated with black hole mass. Even more 
firm confirmation would come from objects with intermediate 
black holes mass, filling the big gap in Fig. [Tfo). One category 
of sources potentially useful for testing this is the ultra-luminous 
X-ray sources (ULX) with possible masses of hundreds of Mq. 
Alas, they are most likely in the very high spectral state, with soft 
X-ray emission dominated by the disc. Besides, there are no re- 
liable high-frequency power spectra available from ULX. An in- 
teresting source for comparison turns out to be a dwarf galaxy 
NGC 4395, with black hole mass estimates between 0.13 and 
3.6xl0 5 M (Kraemer et al. 1999; Ho 2002; Filippenko & Ho 
2003; Vaughan et al. 2005; Greene & Ho 2006). The luminosity is 
low, L/I/Edd < 0.1 (Vaughan et al. 2005). We analyzed its ASCA 
observations of 2000-05-24 and 2000-05-26 and found the excess 
variance for this source and Cm = 3.21*^ x 10~ 5 Hz" 1 . We show 
the range of masses and Cm for NGC 4395 in Fig.|7j6). They are 
in excellent agreement with our mass- variability relation! We can 
also give an independent estimate of the black hole mass in NGC 
4395 based on the best-fitting value of C: M w 3.9 x 10 4 M . 



7 DISCUSSION 

We have analyzed all available hard-state data from seven Galac- 
tic black hole systems and found that the amplitude of the high- 
frequency tail, Cm, is roughly constant for a given source, chang- 
ing by no more than a factor two. There is no apparent depen- 
dence of Cm on luminosity or hardness ratio. This contrasts QPO 
or break frequency behaviour, which typically show a strong cor- 
relation with luminosity. The scale of change of Cm is also much 
less than the observed span of QPO/break frequency, which can be 
one order of magnitude within the hard state. This makes Cm much 
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Figure 7. Dependence of Cm on black hole mass. Data in the upper left panel show the results of this paper for X-ray binaries. Black hole masses, with 
references, are listed in Table[3] The lower right panel shows the sample of Seyfert 1 galxies from N06. Black hole masses are obtained from reverberation 
method (Peterson et al. 2004). The lower left panel show an overview of stellar and supermassive black holes. The red cross in the middle of the diagram 
represent NGC 4395. The solid red diagonal line shows the best-fitting relation Cm — C/M, with C = 1.25 Mq Hz -1 . The green dashed line corresponds 
to the soft state of CygX-1, with C = 2.77 Mq Hz" 1 . 



more invariant feature of a given black hole and a robust estimator 
of its mass. 

There are, certainly, some departures from this rule. Differ- 
ent outbursts of the same transient can have slightly different high- 
frequency tail amplitudes. Most of the sources show an increase in 
Cm below count rate of ~20 s _1 per PCU. This might be attributed 
to systematic instrumental effects at low count rates. Due to large 
errors this has a negligible effect on our results. But, if this effect 
extends to higher count rates, than the difference between 2000 and 
2005 outbursts of XTE Jl 1 18+480 (increase of Cm below 40 s~\ 
see Fig.|2j could be of instrumental origin. 

XTE J1550-564 gives a much clearer exception to this rule. 
The first outburst in 1998 had a very short initial hard state with 
the high-frequency tail amplitude much lower than in the next out- 
bursts. In contrast, the three hard-state only outbursts in 2001, 2002 
and 2003 produced very consistent results. The explanation of this 
phenomenon might lie in the stability of the accretion flow. The 
spectral state during the onset of the 1998 outburst was changing 
faster than in any other outburst analyzed here and took a differ- 
ent track on the colour-colour diagram (Done & Gierlihski 2003). 
It is possible that the accretion flow during such a rapid transition 
was in a different state than in slower transitions, possibly due to 
higher ionization state of the reflector (Wilson & Done 2001; Done 



& Gierlihski 2003). Clearly, this also affected its variability proper- 
ties, decreasing the high-frequency power. This particular hard state 
was different in many aspects, so we rejected it from our sample. 
The more stable hard-state only outbursts give probably a better 
estimate of Cm- Similarly, the persistent source, Cyg X-l, gave a 
very stable and robust high-frequency power. 

Thus, fast transitions, characterized by an unstable and 
quickly changing accretion flow, can apparently break the Cm — 
C /M law. In Table [2] we show the approximate duration of the 
initial hard state after the onset of the outburst. XTE J 1650-500 
showed a rather quick transition, so might have suffered from a sim- 
ilar instability, as XTE J1550-564 in 1998. On the other hand, the 
tail amplitude was very stable and returned to the same level in the 
hard state at the end of the outburst (see Fig.(2]l. The 2004/2005 out- 
burst of GX 339-4 had a much longer hard state than the 2002/2003 
outburst, so perhaps it was a better estimate of the high-frequency 
power. 

In this work we assumed the constant power-law high- 
frequency tail with the spectral index of -2. This is not necessarily 
true and detailed fits to PDS of X-ray binaries, where statistics is 
sufficient, show indices between 1.5 and 2.0. We note that the par- 
ticular choice of the spectral index does not affect the overall result 
of our paper, as different index would only introduced a constant 
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offset in Cm- Scatter in spectral indices from one observation to 
another might introduce some additional uncertainty, though this 
would be very difficult to estimate for AGN. 

The luminosities of the Seyfert 1 sample used in this paper 
are generally low, <5 per cent of LEdd, except for 3C120 and 
NGC 7469, which are significantly brighter. The 3-10 keV photon 
power-law indices are Y ~ 1.5-1.8 (Nandra et al. 1997). This sug- 
gests that these sources are in the hard X-ray spectral state. On the 
other hand, the spectral index of Comptonized component alone is 
not enough to establish the X-ray spectral state. Many BHB show 
soft-state disc-dominated spectra with a flat (F ~ 2) power-law 
Comptonized tail, while measuring the intrinsic spectral index in 
AGN is not straightforward due to presence of complex absorption 
and/or reflection (Gierlihski & Done 2004). We cannot rule out (in 
particular for brighter AGN) that some of the sources in the sample 
are actually in the soft X-ray spectral state. Despite that the mass- 
variability amplitude correlation seems to hold well for all of them. 
Some of the objects are more consistent with the soft-state Cyg 
X-l line in Fig. |7Jc). On the other hand, the brightest 3C120 and 
NGC 7469 lay below the hard-state line. Clearly, the dependence of 
high-frequency amplitude on luminosity and spectral state in AGN 
requires further studies. 

Our analysis of the very high (steep power law) state shows 
that after a simple bandwidth correction the tail power is consis- 
tent with the hard-state results. This is very encouraging, show- 
ing that perhaps the universal shape of the high-frequency tail is 
the inherent property of Comptonization, regardless of the spectral 
state. Narrow-line Seyfert 1 galaxies (NLS1) are most likely the 
supermassive counterparts of Galactic sources in the bright very 
high state (e.g. Pounds, Done & Osborne 1995). N04, and later 
Nikolajuk, Gurynowicz & Czerny (2007) showed that the excess 
variance from the NLS 1 sample is by a factor ~20 larger than ex- 
pected from Cm = C/M correlation established for moderate- 
luminosity AGN. NLS 1 have also systematically higher break fre- 
quencies for a given mass, so the luminosity dependence of the 
break frequency (M c Hardy et al. 2006) can perhaps be related to the 
increase in Cm ■ One explanation could be a bandpass effect. The 
seed photons in X-ray binaries in the very high state are at ~ 1 keV, 
while in AGN they are at much lower energies, < 10 eV. Therefore, 
we see a different part of the Comptonized spectrum in X-rays. 
Gierlihski & Zdziarski (2005) showed that the variability ampli- 
tude strongly increases with energy in the very high state, while it 
is almost constant in the hard state. As we see higher orders of scat- 
tering in NLS 1 than in the very-high-state stellar mass black holes, 
we expect higher variability in the former ones, as observed. An- 
other possible explanation is contribution from complex absorption 
to the observed variability, which can significantly increase the rms 
at energies ~l-2 keV (Markowitz et al. 2003a; Gierlihski & Done 
2006). We expect more complex absorption from bright NLS1 
sources with strong outflows. The low- and moderate-luminosity 
AGN are less likely to be affected by this additional variability. 

Fig. |7] shows that the amplitude of the high-frequency tail in 
the PDS scales very well with the black hole mass. The corre- 
lation holds for over seven orders of magnitude. Thus, the high- 
frequency power in accreting black holes in the hard spectral states 
seems to be universal. This part of the PDS corresponds to the 
shortest timescales producing power in the accretion flow. Fea- 
tures (e.g. QPOs) occasionally observed at higher frequencies (e.g. 
Strohmayer 2001) have much less power. We do not understand the 
origin of rapid X-ray variability from accretion flows very well. It 
might be created by fluctuations propagating in the accretion flow 
(Lyubarskii 1997). The closer to the centre, the shorter characteris- 



tic timescales of fluctuations. Inevitable, there is a final barrier for 
propagation, the last stable orbit, which acts as a low-pass filter, 
removing all frequencies higher than a certain limit. The limiting 
frequency is inversely proportional to the radius of the last stable 
orbit and hence to the black hole mass. With a certain shape of 
the filter any initial spectrum of fluctuations will be truncated to 
a similar shape at higher frequencies. This might (at least qualita- 
tively) explain the observed universal shape of the high-frequency 
tail (Done, Gierlihski & Kubota 2007). 

Obviously, the size of the last stable orbit depends not only 
on the black hole mass, but also on its spin. Higher spin would 
shrink the last stable orbit as if the black hole mass was smaller. 
This would increase Cm by factor 4.8 for a maximally spinning 
Kerr black hole with respect to a Schwarzschild one. Alas, the 
black hole spin is notoriously difficult to measure (compare, e.g., 
McClintock et al. 2006 and Middleton, Done & Gierlihski 2006). 
The potentially highly spinning GRS 1915+105 has never been ob- 
served in the low-luminosity hard spectral state. Another culprit is 
XTE J 1650-500, in which a broad iron line, suggesting high black 
hole spin, has been reported (Miller et al. 2002, but see also Done 
& Gierlihski 2006). However, as one can see in Fig.[7]CM of this 
source is well below the C /M line, so its high spin does not seem 
to be supported by our data. On the other hand the mass of black 
hole in this source is not very well established, so we cannot make 
any strong statements about it. Fig. [7] shows that the scatter in Cm 
for all sources is only factor two, so we do not expect large scatter 
in black hole spins in the sample. 

Yet another factor that might influence the results is the in- 
clination of the disc with respect to the observer. If rapid X-ray 
variability is produced in flares or fluctuations corotating with the 
disc, then it is affected by Doppler effects for highly inclined discs. 
Zycki & Niedzwiecki (2005) calculated these effects and predicted 
that high inclinations would give rise to strong increase in the high- 
frequency power. However, the additional signal appears above 
~ 100 Hz and would not be easily detected in PDS. We calculated 
Cm from 10-128 Hz band, so our results should not be affected by 
the inclination. 



8 COMPARISON WITH BREAK-FREQUENCY SCALING 

An alternative approach of linking variability with black hole mass 
is the mass-luminosity-break frequency scaling, in, — AL^ ol /M c 
(M c Hardy et al. 2006; Kording et al. 2007), where the three observ- 
ables: mass, luminosity and break frequency form a 'fundamental 
plane' along which the BHB and AGN data correlate. 

We would like to point out few advantages of the method pro- 
posed in this paper over the break-frequency scaling. Generally, it 
is easier to find variability amplitude than the break frequency, as 
the latter requires modelling of the power spectrum (but see Pes- 
sah 2007). The amplitude method does not involve luminosity or 
accretion rate, so it doesn't require distance to the source, which in 
case of many Galactic black holes is poorly constrained. Another 
advantage of this approach is its simplicity, as it relates amplitude 
of variability with black hole mass directly, with just one scaling 
constant. The break-frequency method requires three independent 
constants. This suggests that high-frequency variability amplitude 
is more fundamental in nature. 

One of the drawbacks of the amplitude method is its limitation 
to the hard X-ray spectral state in BHB. Another state dominated 
by Comptonization, the very high (steep power law) state is po- 
tentially useful, but its application to bright AGN requires better 
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understanding of energy dependence of rms. Soft-state spectra of 
BHB are strongly diluted by the (stable) disc and good-statistics 
high-energy data, required to establish Cm reliably, is not avail- 
able (except for Cyg X-l and, perhaps, GRS 1915+105). But most 
of the AGN spectra in the 2-10 keV band are dominated by Comp- 
tonization, so this method might be valid in the soft state. The same 
problem seems to affect the break-frequency scaling. Kording et al. 
(2007) point out that their method is mostly limited to the hard 
state, as measuring and defining the break frequency in soft and 
very high states is very difficult. 

Both break-frequency and amplitude correlations require a 
shift in the relation when applied to Cyg X-l in the soft spectral 
state (Kording et al. 2007), though this is rather difficult to extend 
to soft states of other BHB. 

The high-frequency amplitude depends, to some extend, on 
energy, while the break-frequency is energy-independent. Another 
potential disadvantage of the amplitude method is additional vari- 
ability introduced by ionized smeared absorption or reflection 
in some sources, which is most pronounced around 1-2 keV 
(Markowitz et al. 2003b; Gierlihski & Done 2006; Crummy et al. 
2006). 



9 CONCLUSIONS 

Black hole X-ray binaries have a universal shape of the high- 
frequency tail (above the break frequency) in their PDS, as il- 
lustrated in Fig. [T] Though the exact shape of the tail is not 
easy to establish, it can be approximated by a power law, P v = 
Cm{v /vo)~ 2 ■ The amplitude, Cm of the tail is remarkably con- 
stant for any given BHB in the hard state, regardless of the lumi- 
nosity. When extended to supermassive black holes in moderate lu- 
minosity AGN, the tail amplitude scales very well with black hole 
mass, Cm = C/Al. The best-fitting value of the scaling constant 
from our sample of BHB and AGN is C = 1.25±0.06 M Q Hz" 1 . 
This method can be applied to estimate black hole masses in many 
AGN. If the universal shape of the high-frequency tail is an inher- 
ent property of Comptonization, it might be applied to other spec- 
tral states. We speculate that the constancy of the tail is an imprint 
of the last stable orbit around the black hole. 
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